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ABSTRACT 

We present multicomponent solar wind models self-consistently incorporating the contribution from 
dissipationless, monochromatic, hydromagnetic (with angular frequencies w well below ion gyro- 
frequencies), toroidal Alfven waves, which are coupled to the flow through the wave-induced pon- 
deromotive forces. Protons and alpha particles are treated on an equal footing, and the wavelength 
is not assumed small compared with the spatial scales at which the solar wind parameters vary. We 
find that the non-WKB effects are significant, for the fast and slow solar wind solutions alike. In 
comparison with their non-WKB counterparts the WKB ones are more effective in accelerating the 
solar wind inside the Alfven point, producing significantly enhanced ion fluxes and considerably re- 
duced alpha abundance in the inner corona. Only when to > 3.5 x fO~'^ (1.5 x 10~^) rad s~^ can the 
waves in the fast (slow) winds be adequately described by the WKB limit. Moreover, while the Alfven 
waves tend to reduce the magnitude of the proton-alpha speed difference \Uap\ in general, different 
mechanisms operate in two different regimes separated by an ujc ^ several x 10~^ rad s~^. This lOc-, 
defined by Equation (|16p . is closely related to the time required by a solar wind parcel to traverse 
an Alfven radius with the speed of center of mass evaluated at the Alfven point. When u > lUc, the 
fluctuations are wave-like and tend to accelerate both ion species, thereby losing most of their energy 
by doing work on ion flows; whereas when to < Uc, a quasi-static behavior results: the fluctuations 
may act to accelerate the slower flowing ion species but decelerate the faster moving one in a large 
portion of the computational domain, and only a minor fraction of the wave energy flux injected at 
the base is lost. The fluctuations with the lowest frequency are no less effective in reducing \Uap\ than 
the WKB waves: in the slow solar wind solutions, they may be able to quench a significant \Uap\ with 
base amplitudes as small as 4 km s^^. The consequences of ojc on the velocity fluctuation spectra of 
protons and alpha particles, which are likely to be obtained by future missions like Solar Orbiter and 
Solar Probe, are discussed. 
Subject headings: waves — Sun: magnetic fields-solar wind-Stars: winds, outflows 
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1. INTRODUCTION 

Alfven waves have both observational and theoretical 
consequences in solar wind studies. A salient feature 
in the measured solar wind fluctua tions in interplane- 
tary space (see the reviews by, e.g..lTu fc MarschI 119951 : 
iGoldstein et "all [19951 lB"runo fc Carbo^ [2f)n5h . Alfven 
waves may as well account for the Faraday ro tation 
measurements inside 10 Rq (jHollweg et al.l 11982) and 
the non-thermal broadening of a numb er of ultravio- 
let lines measured belo w ~ 5 Rq (e.g., iBaneriee et al.l 
ll998HEsser et al.lll999[ ). It is noteworthy that although 
the hourly-scale fluctuations seem to be more frequently 
studied, the fluctuation spectrum measured by Helios 
nevertheless s pans a broad frequ ency range from 10^^ 
to 10^2 s-i (jTufc MarschI I1995D . On the theoretical 
side, it was recognized even before their identification 
that Alfven waves may provide a ponderom otive force 
that accelerates the solar wind (|Parkea[l963). Further- 
more, Alfven waves may be damped, their energy being 
converted into both thermal and kinetic energies of the 
flow. One such damping mechanism that received much 
attention is based on the idea that a turbulent cascade 
towards high parallel wavenumbers transfers the wave 
energy from the low-frequency hydromagnetic regime to 
the ion cyclotron one, where the energy can be readily 
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picked up by ions t hrough cyclotron resonan ce (see the 
extensive review by iHollweg fc Isenberg||2002l) . 

Most solar wind models that incorporate the contri- 
bution from Alfven waves have been formulated in the 
WKB limit, where the wavelength is assumed much 
shorter than the spatial scale at which the flow param- 
eters vary. On the one hand, this considerably simpli- 
fies the mathematical treatment. On the other hand, 
some measured features of the solar wind fluctuations 
do posses a number of WKB characteristics. For in- 
stance, the Helios measurements revealed that during 
hourly-scale Alfvenic activities in both the fast and slow 
wind, protons are heavily perturbed, where as the velocity 
of alpha particles suffers little distortion (|Marsch et al.l 
119811: iMarsch et al.|[l982allbl L This discrepancy in the re- 
sponses of different ion species to Alfven waves is consis- 
tent with the WKB theory considering that the proton- 
alpha-particle differential speed is of the order of the lo- 
cal Alfven speed. In addition, for the extensively studied 
hourly-scale fluctuations in the fast solar wind around 
1 AU, there tends to exist a high correlation between 
the perturbed magneti c field 5h and plasma velocity ^v 
(jBelcher fc Daviall971|) . As a matter of fact, if examin- 
ing the nor malized cross- hehcity spectrum such as given 
by Fig. 2-4 in iTu fc MarschI (|1995"), one may find that the 
high 6v — (5b correlation is present across the whole iner- 
tial range for fluctuations in the fast solar wind through- 
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out the inner heliosphere explored by Helios. 

As has long been recognized (iHeinemann &: Olberlj 
Il980t) (see also [Goldstein etaLl ([19861)), the apphcabil- 
ity of the WKB approximation is questionable near the 
coronal base where the solar wind is highly inhomoge- 
neous and the phase speed of Alfven waves is high. There 
have been several different approaches to assess the finite 
wavelength effect on Alfven waves. For instance. iHoUwea 
([l973') extended the WKB analysis to higher order and 
found that the finite-wavelength correction may lead to 
a reduction of the wave force below 10 Rq for typical so- 
lar wind parameters. Anoth er approach was adopted by 
IHeinemann fc OlbertI (|1980l ) who developed a formalism 
that is valid for small-amplitude Alfven waves with arbi- 
trary frequencies. This formalism has the advantage that 
outward and inward propagating waves can be explicitly 
separated, the coupling be tween the two can be inte r- 
preted as reflect ion (see alsolHollw eg fc IsenbersI (12007])). 
It was shown bv lCranmer fc Vciri B allcgooiicn (2005) that 
this formal ism is essentially equivalent to the form de- 
veloped by IVellil (jl993^ who expressed the wave trans- 
port in terms of the Elsas ser variables, and equivale nt 
to the equations derived by IZhou fc Matthaeui ()199Clf ) if 
neglecting the nonlinear terms therein. Note that these 
nonlinear terms, which are essential for the development 
of any turbulence, involve the interaction between out- 
ward and inward propagating waves. The adoption of 
Elsasser variables therefore finds application in the one- 
point closure phenomenology that expresses the damp- 
ing of the solar wind turbulence, which at least in the 
case of solar corona is expected to cascade the wave 
energy towards high perpendicular wavenumbers (e.g. , 
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cent advance following this direction is the construction 
of solar wind models that incorporate the contribution 
of Alfven waves through bot h the wave pressure gr adient 
force and turbulent heating (jCranmer et al.ll2"007D . 

Nearly all the non-WKB analyses were conducted in 
the framework of the single-fluid Magnetohydrodynam- 
ics (MHD). In this case, the need to apportion the ac- 
celeration and heating among different species is cir- 
cumvented. However, for the solar wind, alpha parti- 
cles cannot be deemed minor given their non-negligible 
abundance and the fact that they may flow consider- 
ably faster than protons. In the fast solar wind, the 
proton-alpha speed difference Uap = Ua — Up may be 
up to ~ 150 km s~^ a t 0.3 AU where the proton speed 
Up is - 600 km s'^ ()Marsch et al.lll982aD . Moreover, 
a substantial Uap ~ 100 km s"-'^ was also found by He- 
lios in the sl ow solar wind in th e presence of Alfvenic 
fluctuations (jMarsch et al.lfl98lh . Consequently, alpha 
particles can play an important role in determining the 
evolution of Alfven waves. So far the only study that in- 
cludes alpha particles in a non-W KB analysis for Alfven 
waves seems to be that by iLi fc L i (2007) (hereafter pa- 
per I), who demonstrated that substantial deviations ex- 
ist between the WKB expectation and computed wave 
properties even for the relatively high angular frequency 
Lu = 10""^ rad s~^. For the purpose of developing the for- 
malism, paper I examined only the wave propagation in 
a prescribed low-latitude fast solar wind. However, the 
computations indicate that the wave acceleration may 



alter the proton and alpha particle speeds with a wave 
amplitude as low as 10 km s~^ imposed at the coro- 
nal base. This is particularly true for the lowest fre- 
quency considered, to = 10~^ rad s^^, in which case the 
waves may effectively suppress the proton-alpha differ- 
ential speed in interplanetary space. Therefore, a natu- 
ral extension to paper I is to construct solar wind mod- 
els that self-consistently incorporate the acceleration ex- 
erted on ion species by finite wavelength Alfven waves. 
The aim of this paper is to present numerical models thus 
constructed. 

The ponderomotive forces by Alfven waves may also 
play a part in determining the Helium abundance Uap — 
Ha/ Up in the inner corona through their influence on the 
ion flow speeds. It is well established by in situ measure- 
ments for regions ?" ^ 0.3 AU that the fast solar wind 
corresponds to a hardly varying nap ^ 4.5%, and that in 
the slow win d Ugp is highly variab le, ranging from 0.4% 
to 10% (e.g.. lMcComas et al.ll200Cl[ ). However, Uap in the 
corona is subject to some controversy since it has to be 
indirectly inferred. Concerning the remote-sensing mea- 
surements, while evidence exists that Uap niay not be 
enhanced in corona l holes relative to its value for typi- 
cal fast solar winds (jLaming fc Feldnii an 2003) , there also 
exists evidence suggesting the opposite ( Delaboudinierd 
1999). On the other hand, an enhancement in Uap 
at altitudes of one or several tenths of solar radius 
above the limb is nearly inherent t o all multi-fluid mod- 
els, whether an empirica l heating (jHansteen et al.lll997l : 
iLie-Svcnds en et al.ll2003f ) or the ion-cyclotron resonance 
(JLi 2003) is applied to generate the solutions. In some 
cases the corona may even b ecome Helium dominated 
(e.g.. lLie-Svendsen et al.l[2b03f ). Note that in some mod- 
els, the bottom boundary is chosen to be in the mid- 
chromosphere. As a result, the set of governing equations 
has to be extended to include the neutral components of 
Hydrogen and Helium, and the fluxes of both elements 
in the solar wind show a complex dependence on the 
processes that happen in the chromosphere and transi- 
tion region. Neither these processes nor the ion-cyclotron 
resonance will be considered in the present model, how- 
ever. Instead, to isolate the non-WKB effects on Uap, we 
will make our model as simple as possible by placing the 
bottom boundary at the coronal base, and heating the 
ion species by an empirical energy flux. 

The paper is organized as follows. In section |2l we 
give a brief overview of the equations appropriate for a 
system consisting of non-WKB Alfven waves and a multi- 
component solar wind flow. Some further details on the 
numerical implementation of the model are given in sec- 
tion [3l The numerical solutions are then presented in 
sectional Finally, section [5] summarizes the results, end- 
ing with some concluding remarks. 

2. PHYSICAL MODEL 

The solar wind model consists of three species, elec- 
trons (e), protons (p), and alpha particles (a). Each 
species s (s = e,p, a) is characterized by its density n^, 
velocity Vs, mass m,s, electric charge Cg, and tempera- 
ture Tg. The electric charge is also measured in units of 
electron charge e, i.e., e^ = ZsB with Zg = —1 hy defi- 
nition. Similarly, the ion mass number Ak follows from 
the relation ruk — A^mp (k = p,a). The mass density of 
species s is ps = risms, and the species partial pressure is 
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Ps — UsksTs, where ks is the Boltzmann constant. Since 
the wave frequencies of interest are well below the elec- 
tron plasma frequency, the condition of quasi-neutrality 
is guaranteed, He = np+ZaUa- Moreover, quasi-zero cur- 
rent is assumed, Ve = {np'Vp + Zana'Va)/ne, except when 
the meridional ion momentum equations are derived. 

2.1. General Description 

Several simplifying assumptions were made in paper I 
to facilitate the derivation of the governing equations to 
be used in this study. For instance, symmetry about the 
magnetic axis is assumed throughout, i.e., d/d(j) = in 
the spherical coordinate system (r,9,(f>). It is assumed 
that each species considered may be described by the 
standard five-moment equations and that the electron 
inertia may be neglected {rUe = 0). The time- varying 
multicomponent solar wind is then seen as a system 
where the finite wavelength Alfven waves are superim- 
posed on an otherwise time-independent flow, with the 
waves and background flow interacting with each other 
solely through the wave-induced ponderomotive forces. 
The effect of solar rotation is neglected. As a result, the 
unperturbed magnetic field Bp has no (j) components. 
Let e; denote the unit vector along Bp. A flux-tube co- 
ordinate system can then be defined by the base vectors 
(e/. Gat, e^) where e^r = e^ x e/ (see FigH]). Only purely 
toroidal Alfven waves are considered, i.e., the magnetic 
field B and species velocities v^. (s = e,p, a) may be 
written as 



B = Biei + beA,, v^ = Usei + Wgepf 



u..ej 



(1) 



where the lower-case symbols represent perturbations. 

As detailed in paper I, the mathematical manipula- 
tion centers on the fact that for hydromagnetic waves, 
uj is many orders of magnitude lower than the ion gyro- 
frequency flk = {ZkeBi)/{mkc) {k = p,a), where c de- 
notes the speed of light. As a consequence, the ion veloc- 
ity difference vector has to be aligned with the instanta- 
neous magnetic field. Combined with the assumption of 
quasi-zero current, this leads to 



Mfc = Ue + 75" {Uk - Ue) 



(2) 



where k — p,a. Moreover, the velocity perturbations 
Ws [s — e,p,a) are so small that the mass and energy 
exchange between adjacent magnetic flux tubes may be 
safely neglected. One needs to retain Ws only in the sim- 
plified ion momentum equations to ensure the conserva- 
tion of total momentum. In other words, the system of 
vector equations is allowed to be expressed as a force 
balance condition across the meridional magnetic field, 
together with a set of transport equations along it. For 
simplicity, in the present study the force balance con- 
dition is replaced by prescribing a meridional magnetic 
field to be detailed in section 13.11 and only the equa- 
tions governing the time-independent multi-component 
flow and the transport of Alfven waves are considered. 
In essence, this is a 1-dimensional model in that there 
is only one independent spatial variable, I, the arclength 
along a given magnetic line of force measured from its 
footpoint at the solar surface. 

2.2. Equations Governing the Wave Transport 



The equations governing the transport of Alfven waves 
in a multi-component solar wind are the components 
of the magnetic induction law and the total momentum 
conservation (Eqs.(12) and (16) in paper I). Using the 
total momentum as opposed to individual momenta of 
ion species is necessary to eliminate the terms associ- 
ated with ion gyro-frequencies. The analysis of the two 
equations may then be considerably simplified with the 
introduction of the Fourier amplitudes at a given angular 
frequency w, i.e.. 



[bil,t),usil,t)] = [fe(0,Ms(0]exp(-iwi), 



(3) 



where s — e,p, a. Further introducing two dimensionless 
variables 



^^b/Bi, ri^u^/UA, 



(4) 



where Ua = Bi/^Anp is the Alfven speed determined by 
the total mass density p ~ pp + pa, paper I yields that 
the perturbations obey 



(M|- 






Fi2V, 

F22'n- 



(5) 
(6) 



By doing so the original partial differential equations are 
reduced to two ordinary ones, which involve only the di- 
rectional differentiation along the poloidal magnetic field 
d/dl, denoted by the prime ' for brevity. In addition, the 
combined meridional Alfvenic Mach number Mt is de- 
fined by 



M| = 47r(ppC/2 + p,C/2)/i32 



(7) 



The expressions for the complex-valued coefficients 
^11,^12,-^21 and F22, given by Equations (40a) to (40d) 
in paper I, involve the flow parameters pk and Uk {k = 
p, a), the meridional magnetic field Bi, as well as w. Ob- 
viously, Equations ^ and ^ possess an apparent sin- 
gularity at the Alfven point where Mt = 1, which lies 
between 1 Rq and 1 AU for typical solar winds. 

Given a line of force and the background flow parame- 
ters, equations ([5]) and ([6]) may be solved for the Fourier 
amplitudes for the electron velocity and the magnetic 
field perturbations if u> is also known. The ion velocity 
perturbations are then found from the alignment condi- 
tion ([2]) . By assumption, the feedback from waves to the 
solar wind flow is through the wave induced acceleration 
Ofc (fc = p, a) given by 



Ok ^ (ul) {In R}'- 
(bXk) 



Zk 



db 



ATTUg^mk \ dl 



b^(\nR)' 



B, 



(8) 



where the angular brackets represent the time-average 
over one wave period, which may be readily evaluated 
using the Fourier amplitudes (cf. Eq.(21) in paper I). In 
addition, R = r sin ^ is a geometrical factor evaluated at 
each point along the meridional magnetic field line (see 
Fig[T]). The variable Xk is defined by Xk — flkiwj — 
'Wk){Zjnj /ue), where j stands for the ion species other 
than k, i.e., j — a for k — p and vice versa. 

Due to the presence of Wk in Xk , one may expect that 
the iV-component of the ion momentum equations has to 
be solved. In fact, there is no need to do so because Xk 
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involves only the difference Wj — Wk and may be deter 
mined immediately after Uk and h are found, 

duk 

~dt 

'db 



Xk. 



Uk 



-^ +Uk (Ini?) 



Annemk ' 



-Br 



dl 



b(lnRy 



(9) 



The algebraic manipulation leading to this rather in- 
teresting behavior was first devised by iMcKenzie et al.1 
(JT979) when deriving the expressions for the force intro- 
duced into the meridional ion momeii tum equation by 
solar rotation. iMcKenzie et al.l ()1979[ ) noted that this 
practice is in effect an expansion in terms of a small pa- 
rameter nQ/Q,k, where J7q is the angular rotation rate 
of the Sun. That the technique also applies when hydro- 
magnetic Alfvcn waves are concerned is not surprising, 
given that the Alfven waves in this study are represented 
by the azimuthal twists and there also exists a small pa- 
rameter Lu/Qk- 

Equations ([5]) and ([6]) are analytically tractable in the 
WKB limit where uj is sufficiently high to ensure that the 
background flow parameters vary little within one wave- 
length. To the lowest order, the WKB approximation 
yields that (6^) satisfies the equation 



/,,\ P{Uph- U,n)U^h 

{b") 7^ ^— — const, 



Bf 



(10) 



where Um = (PpUp + PaUa)/ P is the speed of center of 
mass, while the phase speed Uph is given by the disper- 
sion relation 



B- 



EPfe(^p'^-f^*0' = 7^, 



(11) 



with k ~ p,a. Furthermore, the amplitudes of the 
species velocity and magnetic field fluctuations are re- 
lated by 



{Uph -Us) 



B, 



(12) 



where s = e,p, a. Finally, a compact expression is found 
for the wave acceleration a^. 



flfc = 



'"^^{b') 



2B 



(13) 



2.3. Equations Governing the Time-independent Solar 

Wind 

From Paper I, the equation governing the meridional 
speed of ion species k reads 



UkUL 



Pk 



Ukirik 



ZkP'e 



Uenik 



GMq 



(Inr)' 



n-i 



+ -r^—CQ{Uj-Uk) + ak, 



(14) 



where k — p,a. The gravitational constant is denoted 
by G, and Mq is the mass of the Sun. In addition, cq 
is a coeflicient associated with Coulomb frictions, which 
unsurprisingly also involve the ion species j other than k. 
We take the Coulomb logarithm to be 21 when eval uating 
Co via the expression given in the appendix of .Li fc Lil 
(|2006l) . Other equations concern the mass conservation 



for both ion species k = p,a, and energy transport for all 
the species s = e,p, a. For simpficity, the ion heat fluxes 
are neglected, whereas the Spitzer law for the electron 

heat flux q^ is assumed, qe 



kT!^^T' where k = 7.: 



10-'^ erg K-^/2 cm-i s'^ ()SpitzeiJll963 ). 

If the wave acceleration ak is a known function of I, 
the model equations can be solved for the distributions 
along a given magnetic line of force of the densities Uk 
and meridional speed Uk of ion species (k — p, a) as 
well as the temperatures Ts of all species (s = e,p,a). 
However, ak is not known a priori but rather depends 
on the flow parameters themselves, thereby coupling the 
waves and solar wind flow. Further complication comes 
from the finite wavelength effect, which introduces an 
apparent singularity at the Alfven point where Mt = 
1 for the wave equations ([5]) and ([6]). The numerical 
method employed for solving such an involved system of 
equations is described in the next section. 

3. NUMERICAL IMPLEMENTATION 

It is necessary to solve the governing equations de- 
scribed in section [3] numerically for a quantitative anal- 
ysis to be made. To this end, the meridional magnetic 
field configuration should be prescribed, some external 
heat deposition needs to be applied to ions to generate 
the solar wind solutions, and appropriate boundary con- 
ditions need to be supplemented for both the wave and 
fiow equations. The implementation of these issues and 
a description of the method of solution are given in this 
section. 

3.1. Background Meridional Magnetic Field 

For the meridional magnetic field, we adopt a n ana- 
lytical model given bv iBanaszkiewicz et al.l (|1998D . In a 
spherical coordinate system, their Equation (5b) may be 
rewritten as 



ipir,0) = M 



sin^ e 3Q sin^ 9 



4-5 sin 



+ 



K 



ai 



fli + r\ cos 6*1 
■\/r2 + a\ + 2air\ cos( 



where ij: denotes the magnetic flux function, whose con- 
tours delineate the magnetic lines of force. The poloidal 
magnetic field is then given by 



[Br, Bg 



1 



rsmi 



Idip 
r'de'' 



dtp 
dr 



In the present implementation, the model magnetic field 
consists of the dipole and current-sheet components only. 
A set of parameters M = 2.9687, Q = 0, X = 0.9343 and 
oi = 1.5 are chosen such that the last open magnetic field 
line is anchored at heliocentric colatitude = 50° on the 
Sun, while at the Earth orbit, the meridional magnetic 
field strength Bi is 47 and independen t of colatitude 9, 
consi stent with Ulysses measurements ([Smith fc Baloghl 
fT995h . 

The background magnetic field configuration is de- 
picted in Figure [U where the thick contours labeled F 
and S represent the lines of force along which we solve 
for fast and slow solar wind solutions, respectively. They 
are also distinguished by the colatitudes 9e where they 
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Fig. 1. — Adopted meridional magnetic field configuration in 
the inner corona. Here only a quadrant is shown in which the 
magnetic axis points upward, and the thick countours labeled F 
and S delineate the lines of force along which the fast and slow 
solar wind solutions are solved for, respectively. Also shown is how 
to define the geometrical factor R, and the base vectors ej, ejv and 
60 of the flux tube coordinate system (see section |2]TJ. 

reach the Earth orbit Re = 1 AU in Table [U which sum- 
marizes the parameters employed to generate the solar 
wind solutions. Tube F (S), which intersects the Earth 
orbit at 70° (89°) colatitude, originates from 9 = 38.5° 
(49.4°) at the Sun where the meridional magnetic field 
strength Bi is 5.24 (4.66) G. 

3.2. Ion Heating 

The external energy deposition is assumed to come 
from an ad hoc energy flux that dissipates at a constant 
length Id- The resulting total heating rate is therefore 

Q^FE—^expl-- 

where Fe is the input flux scaled to the Earth orbit Re = 
1 AU, BiE is the meridional magnetic field strength at 
Re- The heating rate Q is then apportioned between 
protons and alpha particles according to 



t/p Pp 



(15) 



where x is a constant that indicates how alpha particles 
are favored when Q is distributed, with x = 1 standing 
for the neutral heating: the total energy that goes to ion 
species k is proportional to its mass density p^ {k = p,a). 
The choice of the heat deposition needs some expla- 
na tion. An exponential fo rm of Q was first suggested 
bv lHolzer fc AxfordI (|1970( l and was later employed in a 
large number of studies. The form of Q adopted here 
is slightly different from the original version to ensure 
that it mimics the dissipation of a flux of non-thermal 
energy. The way the dissipated energy is distributed 
rese mbles the mechanism involv i ng ion-cyclotron wav es 
(cf. Ilsenberg fc Hollwed ([Up); IHu fc Habbal (fl99l)). 
Previous computations invoking such mechanisms indi- 
cate that, in the case of neutral heating, alpha particles 
tend to flow slower than protons (see the dispersionless 



case in Fig.l of iHu fc Habball ()1999[ )). Only when alpha 
particles are energetically favored in the corona can the 
modeled Uap be positive in interplanetary space. This 
happens when x > 1- Listed in Table [l] the heating pa- 
rameters are chosen to yield fast and slow solar wind so- 
lutions with realistic ion mass fluxes and terminal speeds 
in the absence of wave contribution. As a result, the in- 
fluence introduced by the non-WKB Alfven waves can 
be deduced. 

3.3. Boundary Conditions and Method of Solution 

The governing equations consisting of the flow part 
(Eqs.(8) to (10) in paper I) and the wave part (Eq.(38) 
in paper I) are solved on a spatial grid extending from the 
coronal base (1 Rq) out to 100 Rq- Both parts need to 
be supplemented with appropriate boundary conditions. 
For the flow part, at the base 1 Rq the flow speeds Up and 
Ua are determined by mass conservation, whereas the 
ion densities as well as species temperatures are fixed. 
These parameters are given in Table [1] where the ion 
densities are given in terms of the base values of the 
electron density n^^b and the alpha particle abundance 
{na/np)b- At the top boundary (100 Rq), all the flow 
parameters are linearly extrapolated for simplicity. On 
the other hand, for a monochromatic wave with angular 
frequency w, the regularity requirement at the Alfven 
point means that the wave acceleration Ok is determined 
only by one parameter, which is chosen to be the electron 

1 /2 

velocity fluctuation amplitude (^uf) at 1 Rq, or dub 
for brevity. To evaluate the non-WKB effect, we also 
construct models with WKB waves with the same base 
flow conditions and identical base wave amplitudes. 

We adopt an iterative approach to solve simultane- 
ously the nonlinear system of equations. With an ini- 
tial guess for the flow parameters, the wave accelera- 
tion Ok (fc = p, a) may be obtained as follows. When 
a non-WKB Alfven wave with a given lu is considered, 
equations ([5]) and ([6]) are first solved analytically at the 
Alfven point, and then integrated both inward to the 
coronal base and outward to 100 Rq using the stan- 
dard 4-th order Runge-Kutta method. One may find a^ 
through Equation ([8]). When constructing models with 
WKB waves, the amplitude of the magnetic field fluctu- 
ation |6| is first found at the base from relation (|12p . and 
then computed via Equation (|10p for the whole compu- 
tational domain. The wave acceleration Ok is readily ob- 
tained from Equation (fT5|) . The flow equations (Eqs.(8) 
to (10) in paper I) incorporating Ok a re then so l ved b y 
using the numerical scheme devised bv lHu et al.l (|1997f) . 
thereby all the flow parameters are updated. To en- 
sure internal consistency, the two steps are iterated until 
a convergence is met. This iterative approach cl osely 
follows those used by [A lazraki fc Couturier! (|197H ) and 
iMacGregor fc Charbonneau (,199^ . where WKB and 
non-WKB Alfven waves are self-consistently included in 
single-fluid solar wind models, respectively. 

4. NUMERICAL RESULTS 

Having described the solution method, we are now in 
a position to answer the following questions: to what 
extent the wind parameters are affected by the non- 
WKB effect in general? and how the Helium abun- 
dance and the proton-alpha particle differential speed 
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TABLE 1 
PARAMETERS USED TO GENERATE SOLAR WIND SOLUTIONS 





Magnetic Eield 


Heat Deposition 


Base Elow Parameters 


Wind 
Type 


^bI") 


Fe Id X 
(erg cm-2 s"!) [Rq) 


(10'' K) (IQS cm-3) 


East 
slow 


70 
89 


1.7 1.6 1.7 
1.2 1.4 4 


1 3 6% 
1.4 3 6% 



are affected in particular? To address these questions, 
we will only vary those parameters characterizing the 
waves, namely, the base wave amplitude Sub, and an- 
gular frequency w. There has been considerable effort 
made to determine ^, the non-thermal motions that con- 
tribute to the broadening of some UV lines, in the so- 
lar u pper atmosphere by using SUMER and UVCS data 
(e.g-lBaneriee et aUfToMIChae et aIlll998l:|Peta1ll999t 
lEsser et al.lll999[ ). Tvpicallv ^ is found to be ~ 30 km s'^ 
at an altitude of ^ 0.02 Rq above limb. Since this 
value corresponds to the line-of-sight component, the 
overall magnitude of the non-thermal motion may be 
■\/2£ for circularly polarized waves (cf . iBanerjee et al.l 
ll998HCranmer fc van Ballegooiienl20"05[) . Given that the 
waves examined here are linearly polarized, Sub will be 
allowed to have a value of up to 40 km s^^. In this sec- 
tion, we will fist examine the fast and then the slow solar 
wind. 

4.1. Fast Solar Wind Solutions 

Figure [H examines the flow properties of the fast so- 
lar wind solutions corresponding to several lo as plot- 
ted by different line styles indicated in Figure [5^. For 
comparison, also displayed are the solution incorporat- 
ing WKB waves (the dash-dotted curves) and that with- 
out wave contribution (long-dashed curves, only shown 
in the right column). The models involving waves all 
pertain to a wave amplitude of Sub = 28 km s~^. The 
left column gives the radial distributions near the base 
of the flow r < 2 Rq of the flow speeds of (a) protons 
Up and (b) alpha particles Ua, as well as that of (c) the 
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the right column gives the radial profiles for the entire 
computational domain of (d) Up, (e) Ua, and (g) the dif- 
ference between the two Uap = Ua — Up. The asterisks in 
the right column mark the Alfven point, where the com- 
bined meridional Alfvcnic Mach number Mt = 1, with 
Mt defined by Equation ([7]). Moreover, for the solutions 
with waves the diamonds correspond to the location tnl 
where the amplitude of the wave-induced magnetic field 
fluctuations equals the background, i.e., (&^) /Bf = 1. 
Care has to be taken when one examines the segments 
at distances beyond those diamonds since when deriving 
the wave equation, we assume that the time-dependent 
solar wind flow may be decomposed into a steady back- 
ground and Alfven waves. One may expect that when 
the wave amplitudes are significant such a decom posi- 
tion is not permitted. However, as first discussed bv lLoul 
()1993f l. the source terms introduced into the momentum 
and magnetic induction equations by wave-induced fluc- 
tuations decrease sufficiently fast with distance asymp- 
totically. Consequently, the first-order wave amplitudes 
are valid provided that the wave amplitude imposed at 
the coronal base is sufficiently small. Nevertheless, the 
portions where ^6^) /Bf > 1 are indicated for clarity. 





ElG. 2. — Multicomponent solar wind models incorporating 
non-WKB Alfven waves. Numerical results for three angular fre- 
quencies and for the WKB wind model are given in different line 
styles as indicated in (a). All models have the same base veloc- 
ity fluctuation amplitude Sui, = 28 km s~^. Radial profiles are 
shown for (a) and (d) the proton fiow speed Up, (b) and (e) the 
alpha fiow speed Ua, (c) the alpha abundance Ua/rip, and (f) The 
proton-alpha differential speed Uap = Ua — Up. Panels in the left 
column give a close-up of the inner corona r < 2 Rq. In (d), (e) 
and (f), the long-dashed curves correspond to the waveless model, 
the asterisks indicate the location of the Alfvenic point, defined 
by Equation JTJ to be the location where the combined Alfvenic 
Mach number Mt = 1 . The diamonds mark the location where the 
wave-induced magnetic field fluctuation amplitude starts to exceed 
the background. 



and we have excluded all solutions for which rjsfL < ta 
so that integrating the wave equation backwards from r^ 
to the coronal base is formally allowed. 

A close examination of the left column of Figure [5] re- 
veals that the non-WKB models are virtually indistin- 
guishable. Actually they differ only slightly from the 
waveless model, indicating that the wave acceleration 
associated with these frequencies plays an insignificant 
role in the force balance for protons or alpha particles 
in the region considered. On the other hand, an obvi- 
ous difference exists between the WKB and non-WKB 
wind models. For instance. Figures ^B^ and [^b demon- 
strate that the WKB one tends to produce a higher flow 
speed for both protons and alpha particles near the base, 
and hence a higher ion flux as the ion densities are fixed 
at 1 Rq. Consider the proton and alpha fluxes scaled 
to 1 AU, {npUp)E and {naUa)E- In the WKB case 
[npUp)E {{naUa)E) is 3.18 x 10* (1.49 x 10^) cm^ s'^, 
whereas in the model with w = 10~^ rad s^^, (npUp)E 
and {naUa)E are 2.6x10* and 1.14x10^ cm^ s^^, respec- 
tively. Furthermore, as shown in Figure [^t, the WKB 
case corresponds to a reduced alpha abundance Uap for 
r < 1.56 Rq, with the local maximum n^p^M decreasing 
from 0.19 attained at 1.47 Rq in the non-WKB models 
to 0.14 at 1.55 Rq in the WKB one. This change may be 
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readily understood since compared with their non-WKB 
counterparts, the WKB waves produce an enhancement 
of the alpha speed C/q larger than that of the proton 
one Up. As a result, at the base of the wind the WKB 
model yields a smaller Up/Ua and hence a smaller Ua/np 
since the ion flux ratio {uaUa) / (upUp) is a constant for 
steady flows. An interesting aspect of Figure [Dd is that 
the alpha flow speed Ua in the WKB model becomes 
smaller than in the non-WKB ones for r > 1.58 Rq, de- 
spite that the WKB wave acceleration exerted on alpha 
particles aw,a is substantially larger than the non-WKB 
ones throughout the region considered (see FiglSb). This 
is because close to the coronal base, the reduced alpha 
abundance na/up in the WKB model gives rise to a lower 
heating rate for alpha particles (cf. Eq. fTS)) ') and there- 
fore a lower alpha temperature. Hence the lowered alpha 
pressure gradient force results in a Ua profile that is less 
steep at r > 1.42 Rq. 

Now turn to the right column, from which one may 
gain a first impression of the overall infiuence introduced 
by the non-WKB effect on the fiow parameters. Take 
the variation of the location of the Alfven point, ta, for 
example. For uj = 10~^, 10"'' and 10~^ rad s~^, ta is 14, 
13.6, and 13.9 Rq, respectively. These values are close 
to rA = 14.6 Rq derived for the waveless model. On 
the other hand, va is 11.9 Rq in the WKB model, re- 
inforcing the fact that the WKB Alfven waves are more 
effective than the non-WKB ones in accelerating the so- 
lar wind flow in the sub-Alfvenic portion. However, Fig- 
ure [2ji indicates that even in the WKB case, the proton 
speed Up is only 648 km s~^ at ta, which is only slightly 
enhanced compared with the model without waves, for 
which Up = 577 km s~^ at the Alfven point. More signifi- 
cant difference arises in the super- Alfvenic portion, where 
protons receive a continuous acceleration from both the 
WKB and non-WKB Alfven waves. Furthermore, the 
proton flow speed does not show a monotonic depen- 
dence on wave frequency. For instance, at 100 Rq, the 
WKB wind model yields a C/p = 887 km s~^ whereas the 
non-WKB models yield Up = 686, 778 and 694 km s'^ 
for UJ = 10~^, 10"^ and 10"'^ rad s~^ respectively. For 
comparison, the waveless model obtains 602 km s^^ for 
Up at the same location. 

Inspection of Figure [2J3 shows that among the non- 
WKB models considered, the one with u = 10~^ rad s~^ 
shows a fundamental difference from the rest as far as 
the alpha flow speed Ua is concerned: this wave tends 
to decelerate rather than accelerate the alpha particles. 
In fact, for the entire computational domain Ua is even 
smaller than that in the waveless model. In contrast, 
similar to the WKB case, the non-WKB models with 
cj — 10"'' and 10"'^ rad s~^ produce a higher Ua com- 
pared with the model without waves. The disparate ef- 
fects on the ion flows of fluctuations with to = 10~^ and 
10~^ rad s~^ signify a transition around some tOc be- 
low which the waves behave in a quasi-static manner, 
a feat ure e xtens i vely d iscussed bv iHeinemann fc Olbertj 
(|1980( ) and[Lo3 (|1993f ) in the case of single-fluid MHD. 
Given a background flow, uJc in the present case may be 
approximated by 

UJc » U^A/{2rA), (16) 

where UmA is the speed of the center of mass evalu- 
ated at the Alfven point. Note that when the non- 



WKB waves are self-consistently incorporated, the RHS 
of Equation (fT6|) may depend on the wave frequency and 
base amplitude. However, it turns out that the depen- 
dence is rather weak, a value of Wc ~ 3 x 10~^ rad s~^ 
applies for all the fast solar wind solutions considered. 
For those uj > uJc, one may expect that the waves will 
be increasingly WKB-like with increasing uj, whereas for 
UJ < ojc, the effects of quasi-static fluctuations are similar 
to those of solar rotation, the zero-frequency solution to 
the wave equation. 

Whether or not a quasi-static behavior results, the net 
effect of Alfven waves on the speed difference Uap is to 
reduce its magnitude, as evidenced by Figure [2f. For 



rad s ^one finds that U„n — 



ap 



IS 



0, where the corresponding values 
^ rad s~^ are 28.2 and 45.5 km s~^. 



instance, for uj — 10 
2.87 km s-^at 100 R, 
for UJ = lO"'' and 10 
respectively. As for the WKB model, one finds that 
Uap — 14.7 km s^^. All these values are considerably 
smaller than 108 km s~^, found at the same location in 
the waveless model. This is not surprising since it is read- 
ily shown that if the flow speeds are entirely determined 
by the wave forces, then 



{ul - U'p) 
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const, 



(17) 



where the factor e — 1 for WKB waves. On the other 
han d, relati o n ()17l) also holds when uj approaches zero 
(cf. iLi fc Lil ()2006l )). the only difference is that e should 
be chosen to be 2. For a wave with an arbitrary finite 
UJ, it is natural to expect that the effect lies between 
the two extremes. Now that the ratio (6^) /Bf increases 
monotonically with increasing r for all the wave models 
considered, all models should yield a Uap decreasing with 
r if Ua + Up does not show a substantial variation, as 
confirmed by the super- Alfvenic portions of the flow. 

The behavior of the flow parameters is further ex- 
plained by Figure [31 which displays the radial profiles of 
the wave acceleration inflicted on (a) protons ap and (b) 
alpha particles Oq., as well as (c) the wave energy flirx 
density i^u,. Numerical results are shown for the three 
angular frequencies considered, and the WKB model. 
Figures [3^ and [31d show that within ~ 2 Rq the non- 
WKB wave acceleration differs little from one another, 
but is considerably smaller than that in the WKB case. 
Specifically, one finds that at 2 Rq, for w = 10~^, 10~^ 
and lO"-^ rad s~^, ap is 384, 375 and 373 cm s~^, 
respectively, while the WKB model yields a value of 
1904 cm s~^. At larger distances, the difference between 
the WKB and non-WKB models is also prominent. For 



instance, at 100 Rq, one finds that a„ 



58.9 cm s 



for the WKB model, while 30.2,42.7, and 21.7 cm s" 
for 
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10^^, 10^'' and 10^^ rad s~\ respectively. The 
nonmonotonic frequency dependence of ap is the direct 
reason for the similar dependence of the proton speed 
Up. Likewise, the reduction of the alpha speed Ua for 



10- 



rad 



relative to the waveless model results 



from the fact that in a large portion of the computational 
domain, 5.7 < v/Rq < 31.1 to be precise, the wave ex- 
erts a negative Oq on alpha particles. For w = 10"'^ 
and 10~* rad s~^, however, the waves tend to accelerate 
the alpha particles, as in the WKB case. Note that in 
relation p?)) . the wave effect in limiting \Uap\ does not 
rely on the sign of Uap. Hence we may conclude that 
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Fig. 3. — Properties of non-WKB Alfven waves in 
self-consistently constructed multiconiponent solar wind models. 
Shown are the radial distributions of the wave acceleration ex- 
perienced by (a) protons and (b) alpha particles, as well as (c) 
the wave energy flux density for three angular frequencies and 
the WKB wind model as given in different line styles indicated 
in (a). All models have a common base velocity fluctuation ampli- 
tude Sui, = 28 km s~^. 



the reduction of \Uap\ is achieved in different manners 
by wave-hke and quasi-static fluctuations. While wave- 
like fluctuations tend to provide overall acceleration for 
both ion species, the quasi-static ones may accelerate 
the slower flowing species but decelerate the species that 
flows faster. It is noted, however, that the conclusion 
is based on the condition that \Uap\ < Ua, which holds 
for bot h the fast and slow solar wind as measured by 
Helios ([Marsch et al.|[l982af) . and is also true for all the 
wave-based solutions examined here. For \Uap\ > Ua 
kinetic instabilities may arise when the proton beta is 
com parable to unity a nd the physics is more complicated 
(e.g. iGarv et al.|[2000D . If \Uap\ is larger still, then from 
Equation (jlip one may deduce that a low frequency elec- 
tromagnetic Alfven instability will set in (cf. Verheest 
[l977() . Both subjects are beyond the scope of the present 
paper. 

To proceed, we note that the governing equations allow 
an energy conservation law to be derived, from which it 
is readily recognized that at large distances r ^ va, the 
gain in the ion kinetic and gravitational potential energy 
fluxes derives primarily from the energy fluxes associated 
with waves and the adopted ion heating. In other words, 



by noting that at r ^ r^ the empirical energy flux has 
been exhausted, one may deduce the overall contribution 
of Alfven waves between 1 Rq and r from the relation 

y f pkUkBiE \ Ul{r) + VL 

k—p,a 

^W{RQ)-Wir)+FE, 



(18) 

where W = {Bie / Bi)Fw denotes the wave energy 
flux scaled to the Earth orbit. Moreover, Vesc = 
y/2GMQ/RQ is the escape speed. Equation (|18p indi- 
cates that if the ion fluxes in different models are sim- 
ilar, then the asymptotic ion speed Uk is largely de- 
termined by the difference of W between 1 Rq and r, 
which hereafter will be taken as 100 Rq. Moreover, at 
large distances Ua is close to Up in models with waves. 
Therefore, a frequency dependence of AW similar to that 
of Up is expected and indeed reproduced in Figure [3J;. 
For instance, AW = W{Rq) - W{100Rq) is found to 
be 0.31, 0.67 and 0.36 (here and hereafter in units of 
erg cm~^ s~^) for uj = 10~^,10~'' and 10~^ rad s~\ 
respectively. On the other hand, AW — 1.75 is con- 
siderably larger in the WKB model. In fact, this value 
more than compensates for the enhancement of the ion 
speeds in that among the examined ion fluxes the WKB 
model also yields the largest value. An interesting fea- 
ture is that unlike the frequency dependence of AW, at 
1 Rq the wave energy flux density for uj = 10~^ rad s~^ 
is the closest to the corresponding WKB value. Specif- 
ically, when scaled to Re, the WKB model yields that 
WIRq) = 2.02, whereas the values for W{Rq) are 1.73, 
0.82 and 0.43 for lu = 10"^ 10^^ and lO'^ rad s-\ re- 
spectively. From the definition of the wave energy flux 
density (cf. Eq.(20b) in paper I), it is easily understood 
that this feature comes largely from the frequency depen- 
dence of the wave-induced magnetic fluctuation b. For 
u = 10-^ 10-4, 10-3 and the WKB model, |5| is found 
to be 0.28, 0.15, 0.082 and 0.32 G, respectively. Now 
examine the fraction AW/W{Rq). One may see that 
only 18.1% of the injected wave flux is lost in the form 
of the work done on the ion flows for oj = 10^^ rad s^^, 
whereas for the WKB model the corresponding fraction 
is 86.4%. As for the model with lu = 10^"* (10"^) ^a^j g-i^ 
this fraction is 81.4% (83.2%), close to that found in the 
WKB case. Therefore one may envision that lUc also dis- 
tinguishes fluctuations in their efficiency of losing energy 
to the solar wind flow, with the quasi-static ones being 
much less efficient. 

One may have noticed that in the cases with uj = 10^"^ 
and OJ = 10-"* as well as in the WKB model the F^ pro- 
files are rather similar. As has been discussed in paper 
I, in the near-Sun regions the F^ profiles for all the ex- 
amined waves roughly behave like F^ oc Bi, indicating 
that the wave energy is diluted only by the expansion 
of the magnetic fiux tube. On the other hand, at large 
distances r :^ ta, the F^ profiles for lu = 10-'^ and 
lO^"* rad s^^ roughly comply with the WKB expecta- 
tion, i.e., Fw ^ r~^, since in the region considered the 
WKB limit applies in these two cases. When it comes to 
the model where lu — 10-^ rad s-^, one may find that 
asymptotically f^ ~ r^^ approximately. 

Figure m expands the obtained results by displaying the 
dependence on the angular frequency lu of some flow pa- 
rameters for waves with four different base amplitudes 
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Fig. 4. — Frequency dependence of several parameters characterizing the multicomponent solar wind solutions sclf-consistently incor- 
porating non-WKB Alfven waves. The four series of models pertain to different base velocity fluctuation amplitudes as indicated in (a). 
In addition, the dashed lines give the waveless model, i.e., Suf, = 0. (a) The proton flux scaled to the Earth orbit {npUp)E- (b) The ratio 
of the alpha flux to the proton one {riaU a) / (jipUp) . (c) The local maximum in the inner corona of the alpha abundance {na/np)jvf ■ (d) 
The location of the Alfven point ta- (e) The proton-alpha speed difference at the Alfven point, Ucip,A- (f) The difference of W , the wave 
energy flux density scaled to the Earth orbit, between 1 and 100 Rq. (g), (h), and (i) The proton (7p,ioo and alpha Ua^im flow speeds 
as well the proton-alpha speed difference f/ap.ioo at 100 -Rq. For comparison, the parameters derived in the corresponding WKB wind 
models are given by the short horizontal bars in each panel. 



Sui, = 4,16,28 and 40 km s^^. For comparison, the 
waveless model is shown by the dashed lines, while the 
corresponding WKB results are given by the horizontal 
bars at the right of each panel. Plotted in Figure H] are 

(a) the proton flux scaled to the Earth orbit {npUp)E, 

(b) the ion flux ratio {naUa)/{npUp), (c) the local max- 
imum of the alpha abundance represented by Uap.M, (d) 
the location of the Alfven point r^, (e) the proton- alpha 
differential speed Uap^A at the Alfven point, and (f) the 
loss of the wave energy between 1 Rq and 100 Rq, de- 
noted by AW. Also given are the values at 100 Rq 
for (g), (h) and (i) the proton and alpha speeds, C/p,ioo 
and C/q^iooj as well as the proton-alpha speed difference 
C^Qp,ioo- Let us consider Figures [3^ to^Jj first. One can 
see that as far as the the wave effects inside the Alfven 



point are concerned, for w > 3.5 x 10 '^ rad s 



the 



relative difference of the computed parameters with re- 
spect to the WKB values is less than 10%. In other 
words, the WKB limit provides an adequate description 
for these high frequency waves. However, for lower fre- 
quencies significant differences tend to appear between 
the non-WKB and WKB models, with the latter being 
always more effective in accelerating the solar wind as 
indicated by the smaller ta (see FigHJi). On the other 
hand. Figure [Dd indicates that the difference between 
the WKB and non-WKB values for the ion flux ratio is 
rather modest, even for Sub = 40 km s""'^ the relative 
difference is < 10.2%. Given that the proton flux is sub- 
stantially enhanced in the WKB case (see FigU^), this 
modest difference in the ion flux ratio means that relative 
to the non-WKB waves, the WKB ones boost the fiuxes 



of both ion species to nearly equal extent. Note that 
this docs not guarantee that the alpha abundance in the 
inner corona is only slightly modified by the non-WKB 
effect. In fact, the local maximum Uap^M, which varies 
from 0.168 to 0.182 for 4 x 10"^ < tj < 2 x lO'^ rad s'^ 
in the series where 6ui, = 40 km s^^, is up to 77.3% larger 
than Uap^M = 0.102 obtained in the WKB model. As for 
the proton-alpha speed difference, panel (e) shows that 
the waves with Sut ^ 16 km s^^ have appreciable effect 
in limiting its magnitude already inside the Alfven point. 
Moreover, the wave with ui = 8 x lO^"' rad s^-'^ seems to 
be the most inefficient in achieving this effect in each se- 
ries. Now consider Figures|4f to|4j, where the wave effects 
on the asymptotic flow parameters are examined. It can 
seen in Figure Sf that in general the asymptotic proton 
speed Up is correlated with the wave energy loss AW, 
a consequence of relation (fTS]) . However, the correlation 
disappears for the lowest wave amplitude in the segment 
w < 5 X 10~^ rad s~^. This is not surprising since in this 
case the wave energy loss AW is too low to play a role. As 
a matter of fact, in this frequency range AW is < 0.93% 
of Fe, the energy flux associated with the ion heating. 
However, despite this, the waves show appreciable effects 
in linuting the proton-alpha speed difference (see Fig[4j). 
Specifically, in the waveless model, Uap^ioo = 108 km s~^, 
whereas in the model with uj = 10^^ rad s~^, Uap,ioo is 
59.4 km s~^. The reduction is seen to be achieved by a 
reduction in Ua and an increase in Up (Figsll^ and|4ji). 
Comparing Figure [3J with 0^, one can see that between 
rA and 100 Rq, Uap is further reduced. Furthermore, the 
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feature that the waves with lo = 8 x 10^ rad s^ cor- 
respond to the largest Uap for a given series also shows 
up in Figure UJ. At the higher and lower ends of the 
frequency range, Uap^ioo nearly vanishes for wave ampli- 
tudes Suf, > 16 km s~^. 

4.2. Slow Solar Wind Solutions 

Let us now move on to discuss the slow solar wind, 
or rather, how Alfven waves may influence the flow pa- 
rameters in models obtained along tube S indicated in 
Figure [1] using the heating and boundary flow parame- 
ters corresponding to "slow wind" in Table [TJ For the 
chosen parameters, the waveless model yields a proton 
flux of {npUp)E of 3.73 x 10* cm^ s~^ and an ion flux 
ratio of {nalla) / {upUp) = 0.04. A local maximum of 
the alpha abundance Uap.M = 0.185 is found at 1.9 Rq. 
The flow reaches the Alfven point at 17.2 Rq, where 
the proton (alpha) speed is 258 (367) km s^^, the re- 
sulting speed difference being 109 km s~^. At 100 Rq, 
Up^ioo = 290 km s^^and [/„ — 370 km s^^. It is note- 
worthy that the corresponding Uap = 80 km s~^ is 
not unrealistic for slow solar winds, even larger values 
have been found by Hel ios 2 when approaching perihe- 
lion ( "Marschet anll981[) . 

Figure [5] presents, in the same format as Figure [H 
the frequency dependence of several flow parameters ob- 
tained in four series of solar wind models with different 
Sui,. In addition, the WKB and waveless models are also 
given. A comparison with Figure [4] indicates that nearly 
all the features in Figure [5] are reminiscent of those ob- 
tained for fast solar wind solutions. In particular, Fig- 
ures [IK and EJi indicate that the non-WKB waves pro- 
duce a proton flux up to 32.3% lower, a value for the 
location of the Alfven point up to 31.2% higher than 
their WKB counterparts. As in the case of fast wind so- 
lutions, the WKB and non-WKB Alfven waves show only 
a modest difference in the derived ion flux ratio, indicat- 
ing that the larger wave accelerations in the WKB case 
boost the fluxes of both species to similar degrees (see 
FiglSh)- However, for the local maximum alpha abun- 
dance Hap^M in the inner corona the non-WKB values 
may be up to 91.5% larger (FiglSt). Furthermore, Fig- 
ure [5^ indicates that for Sub > 16 km s~^ both WKB 
and non-WKB waves may substantially reduce Uap in- 
side rA- Note that the waves with uj ^ 2 x 10^^ rad s^^ 
are the most inefficient in achieving this effect, and the 
same feature persists in Figure [H, which shows that Uap 
at 100 Rq is no larger than 27.9 km s~^ for all mod- 
els with 6ub > 16 km s""'^. Interestingly, at this dis- 
tance, the waves with the lowest frequency may also be 
effective even with the lowest amplitude. In fact, for 
cj = 8 X 10~^ rad s~^, the waves with Sub — 4 km s~^ 
yield that Uap.ioo = 24.6 km s"""^, even though the corre- 
sponding wave energy loss is merely 0.013 erg cm~^ s~^, 
amounting to 1.1% of the empirical energy flux exhausted 
to heat ion fluids (cf. FiglSf). In this particular case, 
the wave energy flux plays no part in the solar wind ac- 
celeration, rather waves provide a net acceleration for 
the protons but a net deceleration for the alphas (cf. 
FigslS^ and[5h). Moreover, Figures [5^ and[5h indicate 
that for dub = 16 km s~^, compared with the waveless 
values Ua for w = 8 x 10~^ rad s~^ is barely larger, 
but Up shows a significant increase by 95 km s~^. On 



the other hand, for the WKB and the majority of the 
non-WKB waves, both Up and Ua are enhanced rela- 
tive to the waveless values. This once again reinforces 
the notion that different mechanisms operate at differ- 
ent frequency regimes. A further comparison with Fig- 
ures |4] indicates that some quantitative differences arise 
in Figure[5]nonetheless. One may see that the non-WKB 
values approach the corresponding WKB one at a lower 
frequency. To be more specific, the WKB limit provides 
a description with relative differences within 10% for the 
wave behavior as long as a; > 1.5 x lO^'^ rad s^^. This 
happens in conjunction with the lowering of the critical 
value LOc below which the waves behave in a quasi-static 
rather than a wave-like way. For the solutions examined, 
a value of 1.3 x 10~^ rad s~^ may be quoted for ujc- 

5. SUMMARY AND CONCLUDING REMARKS 

This study has been motivated by the observa- 
tional facts that Alfve nic fluctuations exi s t in both fast 
and slow solar winds (jMarsch et al.l Il981t iMarsch et al.l 
|l982a ,b). and for the majority of the measured fluctua- 
tions the corresponding frequency is too lo w to allow the 
short wavelength WKB fimit to apply (e.g.. lTu fc MarschI 
19951). Spec ifically, using the formulation given in 



Li fc Lil (|2007), we have constructed multicomponent so- 



lar wind models that treat protons and alpha particles 
on an equal footing and that self-consistently incorpo- 
rate the contribution from dissipationless, monochro- 
matic, hydromagnetic (with angular frequencies to well 
below ion gyro- frequencies), toroidal Alfven waves. The 
waves, which smoothly pass the Alfven point where the 
combined meridional Alfvenic Mach number Mt equals 
1 (cf. Eq.([7|)), are coupled to the fiow solely through 
the wave-induced ponderomotive forces. Azimuthal sym- 
metry about the magnetic axis is assumed throughout, 
and solar rotation is neglected. However, no assumption 
has been made that the wavelength is small compared 
with the spatial scales at which the solar wind param- 
eters vary. Starting with a waveless model, for the fast 
and slow solar wind alike, we obtained a grid of models 
corresponding to different lo and base wave amplitudes 

Sub, defined as the amplitude of the electron velocity 

1 /? 
perturbation ^Ug) at 1 Rq. The non-WKB effect is 

then examined, in a systematic and quantitative fashion, 

by comparing the models with that incorporating WKB 

waves with the same 5ub, as shown by Figures [H and [51 

The results may be summarized as follows: 

1. The non-WKB effects are significant for the ma- 
jority of the examined solutions. In comparison 
with comparable WKB waves, the non-WKB ones 
are less effective in accelerating the solar wind in- 
side the Alfven point. As a consequence, non-WKB 
Alfven waves may produce a proton flux up to 32% 
lower, an Alfven point 29% further from the Sun 
than their WKB counterparts for fast solar wind 
solutions. Even though the ion flux ratio increases 
by < 10% from non-WKB to WKB models, in- 
dicating that the WKB waves boost the fluxes of 
both ion species to a similar extent, the maximum 
of the alpha abundance, Uap.M, is considerably af- 
fected by the non-WKB effect, with the non-WKB 
values being up to 77% higher than the WKB ones. 
The differences between WKB and non-WKB mod- 
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Fig. 5. — Similar to Figure 141 but for the slow solar wind solutions. 
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els are more prominent for slow solar wind solu- 
tions, for which the non-WKB models produce val- 
ues for Uap.M up to 92% higher. It is also found 
that the influence associated with non-WKB effects 
tends to decrease with increasing frequency: the 
waves in the fast (slow) winds with lo > 3.5 x 10"'^ 
(1.5 X 10"'^) rad s~^ may be described by the WKB 
limit with an accuracy to better than 10%. 

While the Alfven waves tend to reduce the magni- 
tude of the proton-alpha speed difference \Uap\ in 
general, different mechanisms operate in two dif- 
ferent regimes, separated by a critical frequency 
ujc, which in principle may be different from model 
to model especially when the wave base amplitude 
Sub is large. For the range in which 6ub is varied 
in the examined numerical solutions, however, lUc 
hardly varies with Sub and is found to be ~ 3 x 10~^ 
(~ 1.3 X 10^^) rad s~^ for the fast (slow) solar 
wind models. When to > tOc, the fluctuations are 
wave-like and tend to accelerate both ion species, 
thereby losing most of their energy in doing work 
on ion flows. On the other hand, when uj < lOc, a 
quasi-static behavior results: the fluctuations may 
act to accelerate the slower flowing ion species but 
decelerate the faster moving one in a large por- 
tion of the computational domain, and only a mi- 
nor fraction of the wave energy flux injected at the 
base is lost. The fluctuations with the lowest fre- 
quency are no less efficient in reducing \Uap\ than 
the WKB waves: in the slow solar wind solutions, 
they may be able to quench a significant \Uap\ with 
base amplitudes as small as 4 km s^^. 



Before proceeding, we note that to examine the wave 
effects on the solar wind, an alternative approach will 



be to directly solve the multi-fluid equations where the 
Alf ven waves ar e introduced via boundary conditions 
(cf. lOfmanl [20041 ). In our study, monochromatic waves 
are invoked so that a detailed understanding of the fre- 
quency dependence of the sol ar wind prope rties may be 
found. On the other hand, lOfmanl ()2004l ) introduced 
a broad spectrum of Alfven waves into a resistive and 
viscous solar wind plasma. Consequently, a compari- 
son of his model with ours is n ot straightforward. Nev - 
ertheless, in a previous study (jOfman fc Davilal 11998"). 
where monochromatic Alfven waves were used and a set 
of isothermal single-fluid MHD equations were solved, 
the authors concluded that their results agree with the 
studies of a sing le-fluid isothermal wind as conducted by 
iMacGregor fc Charbonneaul (1994) , who used essentially 
the same approach as ours. From this we may conclude 
that the results obtained in this paper may as well be 
supported by direct numerical solutions, the construc- 
tion of which will be presented in a future publication. 

Another way to extend the current model is to use 
a spectrum of Alfven waves instead of monochromatic 
ones. This will also be more consistent with observa- 
tions which show that in reality, the Alfvenic fluctuations 
span a broad spectrum and do no t possess a preferred 
frequency (e.g. iTu fc Marsclilll995f ). With this caveat in 
mind, the presented study nonetheless provides a better 
understanding of the underlying physics. For instance, 
when a spectrum is taken into account, one may expect 
that net effect will also be that the magnitude of the 
proton-alpha speed difference is reduced, and the reduc- 
tion will have some dependence on the spectral slope. In 
addition, the existence of an Uc of the order of several 
times 10~^ rad s~^ may have significant consequences 
on the spectrum of the velocity perturbations of alpha 
particles, denoted by Pa for brevity. In the present 
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study we have concentrated on monochromatic waves for 
the purpose of presenting a systematic study on the ef- 
fects brought forth by the finite wavelength. However, 
fr om the results obtained in the present paper and paper 
I (|Li fc Lill2007l ). it is possible to deduce some properties 
of Pa in the presence of transverse Alfvenic fluctuations. 
From an observational point of view, the spectrum of 
proton velocity fluctuations, Pp, is less extensively stud- 
ied than that of the magnetic field fluctuations Pb , be- 
cause of the considerably lower temporal resolution of 
plas ma instruments relative t o that of magnetometers 
(cf. iPodesta et aIll2006L [200l . Measuring P^ is even 
more difficult given that the alpha particles are more 
tenuous. As a result, such a study on P„ has yet to 
appear but is certainly feasible for plasma instruments 
on board future missions such as Solar Orbiter or Solar 
Probe. This is because these missions may approach the 
Sun as close as 4 i?©, and therefore the alpha particles 
to be sampled are much hotter and have a much larger 
number density than in the near-Earth regions. 

Let us now picture what Pa in fast solar winds may 
look like in the super-Alfvenic regions (beyond say ~ 
40 i?0 where the Alfvenic Mach number Mr > 3), sup- 
posing that the waves are propagating parallel to a ra- 
dial magnetic field, and that Uap — Ua, Xa = 0.2 and 
Xp — 0.8, where Xk — Pk/p {k — p,Qi). For frequencies 



corresponding to w > coc, the fluctuations are genuinely 
wavelike and one expects that a WKB behavior results. 
It then follows from Equations (25) and (41) in paper I 
that the ratio of the amplitude of alpha to proton velocity 
perturbations r/ « {l — Xp — XpXa/'2)/(2 — Xp — XpXa/'2') ~ 
0.11. On the other hand, for to < uJc, & quasi-static 
manner results, and an analogy is readily drawn with 
the problem of angular momentum transport in a multi- 
compo nent sola. r wind . In this case. Equations (2a) and 
(2b) in lLi et aU (|2007f) . which are equivalent to the zero- 
frequency limit treated in paper I, are more convenient 
to use and yield that 77 « (ppUp) / (paUa) ^ 3 by noting 
that the amplitudes are largely determined by the terms 
associated with Uap- Consequently, in the spectrum of 
the velocity fluctuations of alpha particles Pa , one may 
see an apparent spectral break at Wc around which the 
spectrum shows a steep slope provided that the proton 
velocity fluctuation spectrum is somehow smooth here. 
Whether Paioj) behaves like this in reality remains to be 
tested by future in situ measurements. 
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